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This presentation describes the progress toward developing a general 
method for predicting unsteady heat transfer on turbine blades subject to 
blade-passing frequencies and Reynolds numbers relevant to SSME, such as Illus- 
trated In figure 1. The method employs an Invlscld/vlscous Interactive proce- 
dure (fig. 2) which has been tested extensively for steady subsonic and 
transonic external airfoil problems. One such example Is shown In figure 3. 

The agreement with experimental data and with Navler-Stokes calculations yields 
confidence In the method. The present work extends the technique to account 
for wake generated unsteadiness. The flow reversals around the stagnation 
point caused by the nonuniform onset velocity are accounted for by using the 
Characteristic Box scheme developed by Cebecl and Stewartson. The coupling 
between the Invlscld and viscous methods Is achieved by using a special proce- 
dure, which, with a novel Inverse finite-difference boundary-layer method, 
allows the calculations to be performed for a wide range of flow conditions. 
Including separation. 

Preliminary results are presented for the stagnation region of turbine 
blades for both laminar and turbulent flows. A laminar model problem corre- 
sponding to a flow on a circular cylinder which experiences the periodic pass- 
ing of wakes from turbine blades Is presented to demonstrate the ability of the 
method to calculate flow reversals around the stagnation region. Calculations 
cover a Strouhal number range of 0.1 to 0.2 and a leading-edge Reynolds number 
range of 2xl0 4 to 4xl0 5 . 


BASIC EQUATIONS 

The basic equations are presented In figure 4. As an essential prelimi- 
nary to more extensive calculations Involving complete blades, we consider here 
the development of an unsteady boundary-layer method for calculating the flow 
properties near the stagnation region of a blade where the movement of the 
stagnation point with space and time poses problems. We assume that the exter- 
nal velocity distribution Is represented In dimensionless form by a function 
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(eq. (1)) which allows the variation of the stagnation point and the free 
stream velocity. 

For a two-dimensional, incompressible time- dependent laminar flow, the 
boundary layer equations and their boundary conditions are well known, and for 
conditions with no mass transfer and specified wall temperature they can be 
written in the form shown in equations (2) to (6). 

The calculation of upstream boundary conditions in the (t,y) plane at sortie 

x = x 0 when the conditions at a previous time line are known can introduce 

different problems. To illustrate these difficulties for the case of a moving 
stagnation point, let us consider equation (1). Since u e - 0 at the stag 
nation point by definition, its location based on the external streamlines is 
given by £ s = B(t) . 

Figure 5 shows the variation of the stagnation point with time according 
to the equation with B(t) = 1 + c sin ut, c = 1 , <j - ir/4 . We see that the 
stagnation point is at 2 when t = 2 and at 0 when t = fa, etc. 

If £ s were fixed, we could assume that u - 0 at Es = -1 for all time 

and for all y, but this is not the case. 

It is more convenient and useful to express equations (2) to (5) in a form 
more suitable for computation, lo achieve this, we introduce the dimensionless 
variables t, n, w, m, and G together with a dimensionless stream function 
f(S,n,T) and, with 0 - af/3^, write equations (2) to (5) as shown in equa 
tions (7) to (9). We use the Characteristic Box method to avoid the numerical 
problems associated with flow reversals in the stagnation region and to gener- 
ate the initial conditions on the next time "line." This scheme requires that 
equations (7) and (8) are expressed in terms of new coordinates. For this pur- 
pose, we note the definition of local streamlines and let dr = d^/f ' . If 
the distance in this direction is designated by s and the angle that it makes 
with the T-axis by 6, then the transformed momentum and energy equations (7) 
and (8) and their boundary conditions can be written as shown in equations (10) 
to (13). 


MODtl PROBLEM 

To Illustrate the method we introduce a model problem, which corresponds 
to a flow on a circular cylinder of diameter 0, which experiences the periodic 
passing of the wakes from the turbine blades (see fig. 6). For a time period 
tg the cylinder is subjected to a freestream velocity and for t w it is 

Immersed in a superimposed moving wake which has a motion component up. lhe 
cycle repeats Itself with a so-called blade passing frequency F[E T /( tg + t w ) ] 
and is related to the Strouhal number St by FD/U<». We assume that the stag- 
nation region of the cylinder is subjected to a velocity which varies in space 
and time according to equation (1). (The variables are listed in fig. 6.) 

Calculations are made for two values of Strouhal number taken equal to 0.1 
and 0.2 with the dimensionless minimum wake velocity u m = 1/3 in both 
cases. The computed values of wall heat flux G‘(0) show that they are not 
Influenced by the changes in the freestream velocity and are virtually constant 
for the range of E and t values considered with G 1 ( 0 ) ~ 0.50 for 
St = 0.1 and G'(0) ~ 0.51 for St = 0.2. 
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As shown In figure 7, the computed values of wall shear f“(0) for 
St = 0.1 are significantly Influenced by the changes In the freestream veloc- 
ity which causes flow reversals In the velocity profiles around the stagnation 
point based on the vanishing of the external velocity. The movement of the 
stagnation point and the resulting flow reversals Increase with time and with 
space. For example, the calculations for steady state have the stagnation 
point at 5=0, and, as expected, there Is no flow reversal In either side 
of the stagnation point. At t = 0.05, the stagnation point moves to 
5 = 0.15 but the flow reversals In the velocity profiles continue up to and 
Including 5 = 0.85 as can be seen from the results shown In figure 8. 


CONCLUDING REMARKS 

The computed results show that the numerical procedure Is able to obtain 
solutions for a range of blade-passing frequencies of practical relevance. The 
movement of the stagnation point with space and time and the resulting flow 
reversals around the stagnation point cause no computational difficulties and 
the numerical tests show that the accuracy Is better than required for practi- 
cal problems. 


GENERAL 0BJFCT 1VF 

TO DEVELOP A GENERAL METHOD FOR COMPUTING UNSTEADY HEAT TRANSFER ON TURBINES BLADES 
SUBJECT TO NONUNIFORM ONSET VELOCITY 



APPROACH 

INTERACTIVE BOUNDARY-LAYER THEORY BASED ON THE SOLUTION OF INVISCID AND VISCOUS FLOW 
EQUATIONS BY A NOVEL COUPLING PROCEDURE 


Figure 1. - Approach to rotor wake heat-transfer problem. 
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Figure 2. - Interactive scheme. Viscous 
calculations on the airfoil and in the 
wake are performed by a series of 
successive sweeps to establish new 
boundary conditions for the inviscid 
method . 
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Figure 3. Typical results for single 
airfoi Is. 
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GENERAL FORM OF EXTERNAL VELOCITY FIELO: 


5* Ot.t) . AU.T) [£ - B(t) ] 


BOUNDARY-LAYER EQUATIONS FOR UNSTEADY FLOWS: 


3u 

dx 


+ 


dv 

ay 


o 


in 

at 


+ 


U 4 v |a 
3x 9y 



♦ u 


e dx 



(b 


in 

ay 


) 


3T 

at 


+ u 




yu i_ 

Pr 3y 


«« $ 


y-O, u - v - 0. T- T w (x); 
y - 6. u - u e (x,t), T - T e 


b 



. £t_ 

v Pr t 


(b£") ' + £"6 





£ 


. 3fl 
9C 


1_ 

Pr 


(cG 1 )' + G 1 0 ♦ m(l 


G ) f ' 



£ 


. i£ 

H 


n - o. f - f ■ - g - 0: n - n e . £•-%». G - 1 

"CHARACTERISTIC" FORM OF THE EQUATIONS: 

WRITE THE EQUATIONS ALONG THE LOCAL STREAMLINES 
dx = d£/f 1 

(bf)' ♦ £“0 ♦ |H ♦ w |f - X §£- 

(CG 1 )' + G ' 0 + n( 1 - G)£ 1 - X |2 
T1«0, 0 - f ' - G » 0 , T|» TJ , f 1 « V, 0-1 

- / 1 ♦ (f ' ) 2 . 0- tan _1 f • 

Figure 4. - Description of the method. 
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PRESENT CONTRIBUTION ADDRESSES THE NUMERICAL SOLUTION OF THE BOUNDARY-LAYER 
EQUATIONS FOR A GIVEN PRESSURE DISTRIBUTION AND ITS EVALUATION FOR MODEL LAMINAR AND 
TURBULENT FLOWS. 

FOR NONUNIFORM ONSET VELOCITIES. THE STAGNATION POINT LOCATION t, VARIES WITH SPACE 
AND TIME. I.E. FOR AN OSCILLATING AIRFOIL WHOSE u # NEAR THE LEADING EDGE IS 

2 + 2 

u ~ — (1 + c sin ox) 

6 JT77* 



Figure 5. - Variation of stagnation point with time. 
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Figure 6. - Rotor wake model. 
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. - Flow reversal results for rotor wake model. 
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